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ABSTRACT 



A linear structure can be characterized by its compliance matrix 

C, which is 6x6 S3mimetrical and positive definite and which relates a 

force 6-vector F . to a displacement 6-vector D. by the relation CF . = D,. 

j j j j 

The inverse, S, of C, is called the stiffness matrix and satisfies SD. = 

J 

Fy This thesis deals with the problem of finding optimal values of such 
matrices C and S from experimental determinations of a sufficient number 
of vector-pairs (F^,D^) which are presumed to contain random errors. 

J. E. Brock has introduced this problem area, suggested several dif- 
ferent criteria of optimality, and solved some of the corresponding 
specific problems. This thesis completes the solution to a previously 
unsolved specific problem of this group and contributes computationally 
convenient new solutions to another. Moreover, a computer program, 
originally written for the CDC 1604 has been rewritten, in FORTRAN IV 
Language, as two programs for the IBM System 360 computer, and the 
capability has been significantly augmented. 
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SECTION 1 



INTRODUCTION 

1.1 Source and Statement of the Physical Problem. 

The coordinate system used here is a right handed orthogonal triad, 
whose unit vectors are denoted by e^^, e^, and e^; see figure 1. 

A generalized force will be denoted by F, where 




is a column matrix. Braces £ j* indicate that a column matrix is in- 
tended. The first three elements, F^, F 2 , and F^, are force components 
in the e^, e^, and e^ direction respectively, and the last three ele- 
ments, F^, Fg, and F^, are moment components about these axes. 

Also let a generalized deflection be denoted by D, where 

D . (Di. D2_ 1.3, D3, I.,] 

The first three elements, D^, D^, and are displacements in the 
§ 1 , and e^ direction respectively, and the last three elements are 

rotation components about these axes. 

It is now presumed that a force F is applied to a (suitably con- 
strained) elastic structure, the application being at point P, and the 
corresponding deflection D is observed at P; see figure 1. 

It is postulated that a linear relation holds between the elements 
of F and D, i.e., that there exists a matrix C such that 

CF ^ D 

It is further postulated that matrix C is symmetric and positive definite, 
and that it has an inverse matrix S (S - C which, likewise, is s}mimetric 
and positive definite. 
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Matrix C is called the compliance 
(or flexibility) matrix (of the structure) 
and its elements are called compliance or 
flexibility coefficients. 

Matrix S is called the stiffness matrix, 
and its elements are called stiffness in- 
fluence coefficients. 

Note that matrices C and S depend 
not only on the structure, but also on 
the choice of point P. 

There are basically two approaches for obtaining the compliance 
matrix C, or the stiffness matrix S. One employs elastic theory in an 
attempt to come up with the elements of C or S. Even with simple struc- 
tures it may not be easy to write down the elements of these matrices. 

The theory must be stretched beyond certainty and one obtains only an 
approximation to the correct result. 

The second approach, the use of which leads to this thesis problem, 
employs experimental data rather than elastic theory to come out with the 
best C or S. Because of inevitable random errors present in the data, it 
is impossible to obtain the matric C or S, without using some mathemati- 
cal procedure that is able to minimize the influence of these errors, 
and at the same time to assure that matrices C and S are S 3 mimetrical and 
positive definite. Furthermore the data must contain a sufficient num- 
ber of independent pairs of vectors F and D. 

The question of how many pairs of vectors are required is consider- 
ed in section 7. However in practice, it is anticipated that many more 
pairs will be available than the theoretical minimum. 
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1.2 Brief History of this Problem and Related Problems. 

Professor John E. Brock [1],^ of the Mechanical Engineering Depart- 
ment at the Naval Postgraduate School, introduced this problem. He has 
defined different criteria of optimality and, corresponding to each such 
criterion, has established a matrix equation, the solution of which will 
minimize the corresponding errors and permit determining the desired sym- 
metric and positive definite C or S matrix. 

In the next chapter there will be given the mathematical descriptions 
of four criteria of optimality that Professor Brock has introduced in 
reference [1], and the developments leading to the matrix equations that 
which must be solved for an unknown square matrix. 

A description and an explanation of how to use the computer programs 
which has been developed for solving these equations will be given later. 

The data used to test the programs are not from an actual experi- 
ment. An explanation will be given of how sample data, containing random 
errors, are generated for test purposes. 

In this thesis, matrix C has been restricted to be symmetrical and 
positive definite, as is required in the structural problem under con- 
sideration. Two other related situations will be considered in Appendix 
D. 

In general the different criteria of optimality that will be given 
in the next chapter can be applied to all systems which are described by 
a matrix equation of the form CF = D, where C, F, and D, need only to 
be comformable for matrix operations. 



page 



^Numbers in brackets refer to references listed in the Bibliography 
43. 
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SECTION 2 



MATHEMATICAL DEVELOPMENT OF THE PROBLEM 
2.1 Physical Statement of Criterion One. 

It is assumed that there are n pairs of vector observations F and 

D. 



As criterion one, it is desired to determine C so as to minimize 
the sum of the square of all scalar errors, viz., to minimize 

where = 4-th element of vector 

Now let 



" • ' ^2- ^3- V 

’ - I ®1- 52- ®3 

These are ( ^ x t? ) matrices. 

From equation (1) 






+ C5c^-^iU] <2) 



Expanding the first term of equation (2) 
i ^ + ({CP.-U) +... f (3) 



If all the terms in equation (2), are expanded as in equation (3), 

it is easy to see that these terms are nothing other than the diagonal 

T 

terms of the matrix product: (CU-V) (CU-V) , and therefore it is possible 
to write 
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oc, - U j3(co-v) (cu-v)] 



( 4 ) 



where tr stands for the trace of a matrix and superscript T denotes the 
transposed of a matrix. 

2.2 Alternate Criteria. 




FIGURE 2 

As may be seen from Figure 2, in a one dimensional case, criterion 
one requires minimizing the sum of the square of the vertical errors a^. 
It is natural, therefore to consider, as criterion two, minimizing the 
sum of the squares of the horizontal errors, b^. In the six dimensional 
case that this thesis is concerned with, this means minimizing 

= i I as5j-Fj]j^ (5) 

jif 

Equation (5) is in the same form as equation (1), with the roles of 
Fj and interchanged and using the stiffness matrix S instead of its 
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inverse, the compliance matrix C. 

Using the same approach as before this sum can be written in the 
following way: 

(Xx - [( SV-u/(SV- U) j (6) 

It is possible too, to minimize the sum of the areas of the "error 
triangles", whose legs, in a one dimensional case, are a^ and b^; this 
is the same as minimizing 

(Xj - tn [(cu-v/(u-sv)] (’) 

Now, think also, of minimizing a weighted sum of the square of the 

horizontal and vertical errors ( yS| ^ | This criterion 

number four is the same as minimizing 

oCh =ta[AC<^U-v/(CU-v;-+yS5,(sV-u/(5V-U)] (8) 

where % and are scalar quantities which may be specified so as to 

designate the relative importance of the two separate parts of which it 
is composed. 

From equation (8), without loss of generality, it is possible to 
make j = 1 and , and write equation (8) as 

= th. [(cu-v/(cu-v)+^(sv-u/(sv-u)] (9) 

Brock [1] also suggested a fifth criterion of optimality but since 
no progress has been made toward satisfying this criterion, it is not 
included here. 

2.3 Five Theorems on the Trace of a Matrix. 

In deriving the matrix equation corresponding to each criterion of 
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optimality, the following five simple theorems concerning the trace of 
a matrix will be used. The proofs are simple and are omitted. 



THEOREM 1; If the matrix P and the matrix Q are comformable in 
either order, then each product is square and 
tr (PQ) = tr (QP) 

THEOREM 2: If the square matrices P and Q may be added, then 

tr (P) + tr (Q) = tr (P + Q) 

THEOREM 3: If P is a square matrix, then 

tr (P) = tr (P^) 

THEOREM 4; If P is a matrix and x is a scalar, then 

LsxJ ix 

THEOREM 5: If P is a square matrix and x is a scalar, then 




Note that theorem 1 permits one to write: 

tr (ABCD) = tr (BCDA) = tr (CDAB) , and so on, 

2.4 Derivation of the Matrix Equations which assure the Optimization. 

In sections 2.1 and 2.2 four different criteria of optimality have 
been given. This section devotes attention to the derivation of a matrix 
equation corresponding to each criterion. Solving any such equation 
assures achieving the corresponding optimization, 
a) Equation (4) is 

CX. , = Cl [(C 0 -V)\c. U -V)] 

= Cl (utcu-v'^cu-u^cv+vV) 

T 

Note that instead of C , C has been written since C is symmetric. 

r 

Now let X be any given element of C, and let C* denote ? , then 
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(u''c'^cu+u'cc*u-vt*u-u''c\) 



ztn. [(cUU^ + UUt -UV'^-Vu'"j(C*)] 




where 



P = UU^, and Q = UV^ + VU^ 



are obviously symmetric. 



Thus, if C satisfies the matrix equation 



CP + PC = Q 



( 10 ) 




which will be minimum, since oC\ Is the sum of the squares of real 
numbers . 

Because of the simple relationship between oC% and 0(2^ , it is 
evident that a similar equation applies for criterion two. 

This equation is 



The solution will be to determine the matrices C which make ^ 
stationary; among these will be one which minimizes 0(3. 



SR + RS = Q 



( 11 ) 



T 

where Q is as before and R = W . 



b) Equation ( 7 ) is 




Knowing that S = C and both C and S are s)mimetric, after some 



manipulation 
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c <5 = (u^CU“V'^U-U'"v+V^5v) 

As before, let x be any given element of C, and let C* = 

S* 

Since CS = I, then 
CS* + C*S = 0 
and S* = -SOS. 

Thus 

-Ifi- r tL(u’c*U-V‘^SC*SV) 

= fe[(uu'^-svv'^s)(c*')] 

This will vanish and oC^ will be extremal if S satisfies the matrix 
equation 

T T 

UU = SW S 

After post and pre-multiplying by C 

CUU^C = (12) 

c) Equation (9) is 

°<H = tA[(cU-v]*^(cu-v) + ^ (sv-uj"^(sv-u)j 
= Sl[v\-u'^CV-V*^CU +utcu 
ys(v''5sv-usv-v’su +u^u)] 

As before, let x be any given element of C. 

Thus, using notation introduced previously, viz. : 

P = UU^, R = W^, and Q = VU^ + UV^ 

it may be seen that 
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^ \/c*y - y \ *u + o'" c*c u + u'^c + 



^ (v^*sv + v”^ss*v - u^s *v -v‘^s*u)‘] 

= ^ [(-vu"^-uvVcuu'^4- uu'^c )(c^] + 

yS (svv'^+w'^s -VU"^- uv'^)(s^)] 

= fe [(-a + CP+PC)(C*)+y6(sR^RS-Q)(-sc''s;] 

= 5;, [{ (-a+PC+CP)-^S(sR-fRS-Q)s } (C*)] 

which will vanish if 

PC + CP - Q = y3 S(SR + RS - Q)S (13) 

If the right hand side of equation (13) is set equal to zero, the 
left hand side is precisely the same as equation (10) , which corresponds 
to criterion number one. If the left hand side is set equal to zero, 
the right hand side is nothing other than equation (11), which corres- 
ponds to criterion number two. 
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SECTION 3 



SOLUTION OF EQUATIONS FOR CRITERIA 1 AND 2 
Equations (10) and (11) that minimize criteria 1 and 2 respectively 
are of the same form, where the matrices C, S, P, R, and Q, are all sym- 
metric. The solution by Brock [1], given here, is for equation (10), but 
is applicable to equation (11), with the roles of P and R interchanged. 
The solution in equation (10) will be for the unknown matrix C, where 
in equation (11) is for the unknown matrix S. 

In the physical problem of primary interest, matrices C, P, and Q 
are all 6x6, but, for the sake of generality, it will be presumed that 
they are ra x m. 

In equation (10) 

CP + PC = Q 

there are really only N = % m (m + 1) unknown elements of C. These un- 
knowns can be arranged in a column matrix according to the following 
procedure ; 

Use the notation = % (j - 1) j. 

Given i, j, let t + min (i,j). 

i,j max(i,j; 

For example, with i = 3, and j = 5, then 

t3,5 = t^^3 = % ( 4 )( 5 ) + 3 = 13 

and t^^^ = t ^^2 = % (5) (6) + 2 = 17 

Table 1 gives values of t. , for i,j = 1, 2, 3,. ..,10 

^ > J 

Next, a column matrix W of N elements is formed from the N unknown 

elements c. . of matrix C. Element c. . of C becomes the t. ,th element 
ij iJ 

of W; note that t = N. 

m,m 
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1=1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


j=l 


1 


2 


4 


7 


11 


16 


22 


29 


37 


46 


2 


2 


3 


5 


8 


12 


17 


23 


30 


38 


47 


3 


4 


5 


6 


9 


13 


18 


24 


31 


39 


48 


4 


7 


8 


9 


10 


14 


19 


25 


32 


40 


49 


5 


11 


12 


13 


14 


15 


20 


26 


33 


41 


50 


6 


16 


17 


18 


19 


20 


21 


27 


34 


42 


51 


7 


22 


23 


24 


25 


26 


27 


28 


35 


43 


52 


8 


29 


30 


31 


32 


33 


34 


35 


36 


44 


53 


9 


37 


38 


39 


40 


41 


42 


43 


44 


45 


54 


10 


46 


47 


48 


49 


50 


51 


52 


53 


54 


55 



TABLE 1. Values of t 

* » J 

“■l”l "nI'C'u 'hm} 

The column matrix Z of N elements is formed in a similar way from 

the known elements q, , of symmetric matrix Q. 

J J 

^ *3 

where e . / (i + ^ 

^ f J 

and Kronecker's delta. 

Finally the coefficient matrix P is formed as the sum of m (N x N) 

matrices P^, I * l,2,...,m, each of order N x N, which themselves are 

T 

formed from the (m x m) matrix P = UU , according to the following pro- 
cedure: 

Pj, is formed from P by inserting rows and columns of zeros, an oper 

atlon that is denoted by (*)(K)(L) which means insert L rows of zeros 

after the row in the original matrix P, and Insert L columns of 
th 

zeros after the K column in the original matrix P, The P^ is obtained 
from P by the operation: 
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(a) (*)(0)(M^ 

(b) (*)(I)(I-1) 

(c) (*)(I+1)(I) 

(d) (*) (1+2) (I+l) , etc. 

until such insertions fall outside the range of the N x N size. 

Then: 

yr> 

p - X Pj 

J=l 

and solving the system 

PW = Z 

gives the column matrix W, whose elements can then be rearranged to 
produce the desired matrix C. 

Bickley and McNamee [4] have considered a more general equation of 
the form: 

AZ + ZB = F 

and they give four methods of solution, which they have named: The Big 

Matrix, the Irrational Solution, the Semi-Rational Solution, and A Ration 
al Solution. 

The method given by Larsen [5], and the method given by Brock [1] 
and reproduced here, appear to be additional distinct methods; both, 
however, probably are closely related to what Bickley and McNamee refer 
to as the ”big matrix”, in which A =* B. 

As the results obtained with the method given here is quite satis- 
factory, no further study has been done with the methods given in refer- 
ence [4], but it might be of interest to one working with this problem 
to try suitably specializing Bickley and McNamee *s methods. 
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SECTION 4 



SOLUTION OF EQUATION FOR CRITERION 3 

T 

To solve equation (12), postmultiply by UU J then 

CUU^CUU^ = W^UU^ (14) 

CUU^ = (W^UU^)^ (15) 

T T 

In equation (14), UU and W , both are symmetric, but the product 
T T 

W UU is non-symmetric . 

Brock [6], has solved equation (15), using a highly sophisticated 

computer library subroutine, based on a Laguerre-Hessenberg algorithm, 

T T 

to find eigenvalues of W UU , which are then used in a method described 
by Frazer, Duncan, and Collar [3], in order to find CUU*^ = (w'^UU*^)^, 
from which C may easily be obtained. This procedure is available in a 
computer program for the FORTRAN 63 language to be used in the CDC 1604 
computer. 

Owing to the difficulties in converting this program to FORTRAN IV 
language used in the IBM System 360 computer, recently installed at the 
Naval Postgraduate School^ the following two other procedures were develop- 
ed and have successfully been applied to solve equation (12). 

a) Iterative procedure for finding the square root of a non- 
symmetric matrix. 

Pulay [7], has given a method that determines the principal square 
root of a positive definite symmetric matrix using an iteration pro- 
cedure. 

Using his idea with some modification, a method for obtaining the 
square root of a positive definite non-symmetric matrix has been suc- 
cessfully used and applied to several examples. 

In equation (15) 
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T T T i 

CUU = (W UU 

T T 

let W UU = A 

and = D + B (16) 

where D is a diagonal matrix, whose elements are the positive square 

^ T 

roots of the diagonal elements of A, and B is a matrix formed from CUU 
D, and C is an approximation of the matrix C that will minimize the cri- 
terion number 3. 

fsJ 

The matrix C can be obtained from the equation 

CU = V 
CUU^ = VU^ 

and C = VU^(UU^)"^ 

This is one way of finding the matrix C, suggested by analogy to 
Pulay*s method. However, in the computer program used in this thesis 
to solve equation (15), the matrix 2 is taken to be the matrix C that 
satisfies criterion number 2. This is done in order to save computer 
time . 

From equation (16) 

(D + B)^ = A 
DB + BD = A - - B^ 

2 

This equation can not be solved for matrix B due to the term B , 
but following Pulay*s procedure it can be transformed into a simple 
iteration formula 

“k + . A - 

and the (i,j)-th element of B , can be obtained in a simple manner from 
the formula 
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To obtain iteration procedure must be carried out until 

the difference between two successive approximations of become 

less than a specified allowable error, i.e. 



(B, ) . . - (B, A 

^ k^ij ^ k-l^ij 



where in the actual program is 10 and the Oj- j implies comparing 

element by element. 



After matrix B has been obtained 



C = (D + B)(UU^)"^ 

C = (W^UU^)^(UU^)"^ 



The procedure described above leads to output 3A in the actual 

computer program. Output 3B is the solution of basically the same equa- 

T 

tion, but instead of equation (12) being postmultiplied by UU , it is 
premultiplied resulting in 

UU^C = (UU^W^)^ 

T T 

and using the same procedure as before the square root of UU W can be 
obtained, and therefore the unknown matric C* 

Output 3C is simply the arithmetic average of 3A and 3B. 

b) The second procedure that solves equation (12) uses some 
matrix algebra manipulation, and instead of working with a non-sjmimetric 
matrix , it requires taking the square roots of two different S3mimetric 
matrices. This procedure was called to the attention of the writer by 
Professor Brock. 

T T 

Observe that starting with equation (12), CUU C = W and pre- 
T 

multiplying by U and postmultiplying by U, the following matrix equa- 
tion is obtained 

T T T T 

U CUU CU = U W U 



and 
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U^CU = (U^W^U)^ 



T T 

where U W U is a synnnetric matrix, but as U and V have dimensions (6xn) , 
where n is the number of pair of vectors available, this formulation of 
the problem requires taking the square root to an (nxn) symmetrix matrix, 
with the disadvantage of dealing with a big matrix, for example (1000x1000). 
Accordingly, return to equation (12) 



T T 

CUU C = W 

T 

Let UU = A, and without much loss of generality assume that the 
eigenvalues of A are distinct. In this problem A is a (6x6) matrix. 
Then 






where M is the modal matrix of A 



and L is the spectral matrix of A. 



This is true since A is symmetrical, L is simply a diagonal matrix 

of the eigenvalues of A, and M is a square matrix whose columns are the 

T -1 

corresponding eigenvectors of A. M is orthogonal so that M = M 
2 

Since L is a diagonal matrix, it is easy to find its square 
root, L, Theoretically there are a number of different choices of L 
since either the positive or the negative square root of each eigen- 
values can be chosen, but in this application it is convenient and satis- 
factory to use only the positive values. 

A = ML V = ULUi'^ 

= mllV 

= EE^ 



where E = ML is (mxm) , i.e,, 6 x 6 in the applied problem. 
Let B = 

and equation (12) now takes the form 
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CEE^C = B 



E^CEE^CE = E*^BE = D 



r 

where D is an (m;<m) 


symmetric matrix. Thus, 

T % 

E CE = 



Like matrix A, D can be written in terms of its modal and spectral 



matrices 


D = ml^m^ 
= mlm'^ 



Therefore 



Now 


T T 

E CE = mlm 

e”^ = l”^'^ = l"V 


where L ^ is simply 


the diagonal matrix of the reciprocals of the corres 



ponding elements of L, 



Thus 


-1 T -1 T 

C = ML mlm L M 



This procedure has been written in the FORTRAN IV language to be 
used in the IBM System 360 computer, and in the actual program is output 
3D. 
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SECTION 5 



SOLUTION OF EQUATION FOR CRITERION 4 
The equation (13), that minimize the criterion 4, is nonlinear in 
the unknown matrix C (= S ^) , and an iteration method is used for its 
solution. 

Given , choose a convenient f ^ ^ 
consider the iterative scheme indicated by 

C.P + PC. - Q = f S.(S.R + RS. - Q)S, 

1 I 'l l l i 



Now from equation (13), 



(17) 



where 



and 



S. = C.", 
1 1-1 



C = 
o 



r <■ 

rry 



is the matrix that satisfies criterion 1 
C 2 is the matrix that satisfies criterion 2 
Now, having C^, is determined, from which the right hand side 
of equation (17) can be evaluated, and the following equation can be 
obtained: 



C.P + PC. - Q* 
1 1 



(18) 



where 



Q* = Q + ^ 

Equation (18) is of the same form as equation (10), and can be 

solved using the same procedure. 

The C. obtained from equation (18), is used again to obtain S. 

and a new C.,^ is obtained using the iterative procedure until two 
i+i 

successive approximations of C^ become less than a specified allowable 
error, viz. 



^^i+l^jk " ^^i^jk ^ 
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where ^ in the actual program is taken as 10 
comparing element by element. 



, and Ojj^ implies 



-5 



When the convergence is achieved, ^ is increased and the iterative 
procedure is carried out again. This procedure is followed until finally 
convergence is obtained with ^ = Ji 

The reason for using values of ^ initially less than jg is to 
assure and optimize convergence of the entire process. For example, in 
one study where Ji =0.8, it was found convenient to use the sequence 

f = 0.1, 0.2, 0.3,. ..,0.8. 
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SECTION 6 



DIGITAL COMPUTER PROGRAMS 

Two main computer programs have been written in FORTRAN IV Language 
to be used with the IBM System 360 computer. These programs are called 
CECI 1 and CECI 2, They make use of several subroutines and to assure 
accuracy the main programs and the subroutines employ double precision 
arithmetic. 

In program CECI 1, the following subroutines are used: 



ADDEM 


AINVM 


RIGCM 


BIGSM 


BLANCA 


JACVAT 


MAUSS3 


OUTPU 


PRINM 


RECIP 


REDEM 


SQRERT 


SQROOT 


SQUAM 


TRACM 


TRANM 


ULTMA 


ULTMB 


In the program CECI 2, 


the following subroutines are used: 


ADDEM 


AINVM 


BIGCM 


BIGSM 


BLANCA 


MAUSS3 


OUTPUT 


PRINM 


REDEM 


SQUAM 


TRACM 


TRANM 


ULTMA 


ULTMB 



Both programs make use of a number of functions from the monitor 
library tape which are not listed here. 
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Program CECI 1 solves the equations that satisfy criteria numbers 

1, 2, and 3. The program gives one output for each of criteria 1 and 

2, and four outputs for criterion number 3. For each of these outputs 
the program prints the compliance matrix C, and for each such matrix, 
computes and prints out three criteria numbers, as follows: 

CRITERION NUMBER ONE [(CU - V)^(CU - V)] 

CRITERION NUMBER TWO =tn. [ (SV - U)^(SV - V)] 

CRITERION NUMBER THREE =t>L [ (CU - V)^(U - SV) ] 

Program CECI 2 solves the equation that satisfies criterion number 
4, prints the compliance matrix C, and gives four criterion numbers; the 
first three are the same as in CECI 1 and the fourth is given by: 

CRITERION NUMBER FOUR = t>t [ (CU - V)^(CU - V) + J2 (SV-U)^ 

(SV- U)] 

Most of the matrices are regarded as doubly subscripted arrays as 
follows : 

If there are M rows and N columns, the array in the computer memory 
has dimension (M + 1, N); for example: if the matrix has 6 rows and 50 

columns, the dimension of the array is (7 x 50). This extra row permits 
the identification of the matrix for debugging purposes and gives this 
information to those subroutines which use it for processing or monitor- 
ing purposes . 

Let e. . refer to the array and let a. . refer to the matrix. Then 
i,J i,J 

e - is an integer identifying the array and the matrix, e and e- ^ 

1,1 1,21>J 

are integers M and N respectively. For the remainder of the elements, 

e,,- . = a. .. The other matrices used in the main programs that have 

1+1, J i,J 

the dimensions (N,N) , are used as scratch matrices. They are needed 
when calling some subroutines that have dummy subscripted variables, due 
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to the mode of storage (column-wise) used by the IBM System 360 computer. 
The array for an M-element column matrix has dimension (M + 1), where 
the first element is the integer M. 

The programs can be used in two different ways: (1) for classroom or 
more investigation purposes and (2) for use with data from an actual ex- 
periment. This is done by means of the card called MODE, where MODE = 1, 
2, respectively. 

If MODE = 1, the program will read the matrices C and U from the 
data cards, and generate matrix V = CU. For testing purposes random 
errors are introduced in matrices U and V. The explanation of how these 
errors are introduced will be given in the next section. 

If MODE = 2, the programs will read the matrices U and V from the 
data cards. This data would be obtained from an actual experiment. 

The explanation of how to use these programs, and how the data is 
punched in cards will be given in Appendix A. 

Unless other FORMAT statements are called for, the only printed out- 
2 

put is as follows: 

a) In main program CECI 1, with MODE = 1. 

(1) MATRIX P FOLLOWS 

MATRIX NUMBER 1, 6 ROWS, 6 COLUMNS 
Listing of matrix P in FORMAT (IX, 8D12.4) 

(2) MATRIX BE FOLLOWS 

MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix BE in FORMAT (IX, 8D12.4) 



The nomenclature used in the outputs of the programs CECI 1 and 
CECI 2 is given in Appendix C; see pages 76 and 89. 
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(3) MATRIX EB FOLLOWS 

MATRIX NUMBER 3, 6 ROWS, NN COLUMNS 
Listing of Matrix EB in FORMAT (IX, 8D12.4) 

(4) ERROR IN BLANCA = (number given) 

(5) MATRIX U FOLLOWS 

MATRIX NUMBER 4, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (IX, 8D12.4) 

(6) MATRIX V FOLLOWS 

MATRIX NUMBER 5, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (IX, 8D12.4) 

(7) ELAPSED TIME = (time given in sec.) 

(8) OUTPUT NO. 1, MINIMIZES V SQUARE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER 17, 6 ROWS, 6 COLUMNS 

Listing of desired compliance matrix in FORMAT (IX, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQUALS (value given) 

(9) ELASPED TIME = (time given in sec.) 

(10) OUTPUT NO. 2, MINIMIZES U SQUARE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER 33, 6 ROWS, 6 COLUMNS 

Listing of desired compliance matrix in FORMAT (IX, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQUALS (value given) 

(11) ELAPSED TIME = (time given in sec.) 
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(12) OUTPUT NO. 3A, MINIMIZES TRIANGLE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER XXX, 6 ROWS, 6 COLUMNS 

(XXX may have any large integer value, due to the iteration 
procedure used to obtain output 3A) 

Listing of desired compliance matrix in FORMAT (IX, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQUALS (value given) 

(13) ELASPED TIME = (time given in sec.) 

(14) OUTPUT NO. 3B, MINIMIZES TRIANGLE ERROR 
Same as number (12) 

(15) ELAPSED TIME = (time given in sec.) 

(16) OUTPUT NO. 3C, MINIMIZES TRIANGLE ERROR 
Same as number (12) 

(17) OUTPUT NO. 3D, MINIMIZES TRIANGLE ERROR 
Same as number (12) 

(18) ELAPSED TIME = (time given in sec.) 

(19) TOTAL ELAPSED TIME = (time given in sec.) 
b) In main program CECI 1, with MODE = 2. 

(1) MATRIX U FOLLOWS 

MATRIX NUMBER 1, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (IX, 8D12.4) 

(2) MATRIX V FOLLOWS 

MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (IX, 8D12.4) 

(3) ELAPSED TIME = (time given in sec.) 
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(4) OUTPUT NO. 1, MINIMIZES V SQUARE ERROR 
DESIRED MATRIX IS 

MATRIX NUMBER 14, 6 ROWS, 6 COLUMNS 

Listing of desired compliance matrix in FORMAT (IX, 8D12.4) 
CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQUALS (value given) 

From here the printed out is the same as when MODE = 1. 
c) In main program CECI 2, with MODE = 1. 

(1) MATRIX P FOLLOWS 

MATRIX NUMBER 1, 6 ROWS, 6 COLUMNS 
Listing of matrix P in FORMAT (IX, 8D12.4) 

(2) MATRIX BE FOLLOWS 

MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix BE in FORMAT (IX, 8D12.4) 

(3) MATRIX EB FOLLOWS 

MATRIX NUMBER 3, 6 ROWS, NN COLUMNS 
Listing of matrix EB in FORMAT (IX, 8D12.4) 

(4) ERROR IN BLANCA = (value given) 

(5) MATRIX U FOLLOWS 

MATRIX NUMBER 4, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (IX, 8D12.4) 

(6) MATRIX V FOLLOWS 

MATRIX NUMBER 5, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (IX, 8D12.4) 

(7) ELAPSED TIME = (time given in sec.) 

(8) ELAPSED TIME = (time given in sec.) 

(9) NUMBER OF ITERATIONS = (number given) 
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(10) OUTPUT NO. 4, MINIMIZES V + BETA x U SQUARE ERRORS 
DESIRED MATRIX IS 



MATRIX NUMBER XXX, 6 ROWS, 6 COLUMNS 

Listing of desired compliance matrix in FORMAT (IX, 8D12.4) 
BETA = (number given) 

CRITERION NUMBER 1 EQUALS (value given) 

CRITERION NUMBER 2 EQUALS (value given) 

CRITERION NUMBER 3 EQAULS (value given) 

CRITERION NUMBER 4 EQUALS (value given) 

(11) ELAPSED TIME = (time given in sec.) 

(12) TOTAL ELAPSED TIME = (time given in sec.) 
d) In main program CECI 2, with MODE = 2. 

(1) MATRIX U FOLLOWS 

MATRIX NUMBER 1, 6 ROWS, NN COLUMNS 
Listing of matrix U in FORMAT (IX, 8D12.4) 



(2) MATRIX V FOLLOWS 

MATRIX NUMBER 2, 6 ROWS, NN COLUMNS 
Listing of matrix V in FORMAT (IX, 8D12.4) 

From here the printed out is the same as when MODE = 1. 

Actual printouts for both programs, when MODE = 1, are shown in 
Appendix C. 

The listing of the main programs and subroutines will be given in 
Appendices A and B respectively. 
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SECTION 7 



TESTING OF THEORY AND DIGITAL COMPUTER PROGRAMS 



In this investigation the digital computer programs written to solve 
the matrix equations that satisfy the four criteria of optimality have 
not been tested using actual (mechanical) experimental data. The data 
have been generated by assuming a symmetric positive definite matrix C, 
and a set of vector forces ^F.> in the form of matrix U. From these two 
matrices the corresponding set of vector deflections { ® the form 
of matrix V were obtained as V = CU. 

To simulate experimental data, matrices U* and V* were formed in the 
following way:^ 



U* = U + E(U) 
V* = V + E(V) 



where E(U) and E(V) are matrices of errors. 

Then matrices U* and V* were employed as if they were composed of 
experimental data, to be used in the programs to recover an approximation 
matrix C* for C. 

To form error matrices E(U), and E(V), a mathematical procedure was 
devised to perform some alterations in each element of the matrices, U 
and V, in a random manner with an upper limit to the absolute value of 
the percentage of error being introduced. This procedure has not been 
tested for randomness, but the method used to obtain the alterations 
suggests that the values of the elements in matrices E(U) and E(V) are 
indeed pretty good random errors. 



Note that the use of the asterisk,*, is different in this section 
from its use in section 2. 
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The random errors are generated by a subroutine called BLANCA. 

Because the first row of most of the matrices are used for identifying 
purposes, the subroutine makes use of the element(l, 1) , which is an 
identifying integer called NUMBER. 

Subroutine BLANCA modifies a given matrix A = [ a^j ] by introduc- 
ing random errors E = resulting in a modified matrix A* = A+E. 

One specifies the maximum percentage error to be applied. 

First a number ^ = 3.6327 + (matrix identification NUMBER) is 
formed. ^ is increased for each element by adding a real constant (3,6327) . 
It is not allowed to be greater than 9.999; this is done by subtracting a 
real constant (8.4153) each time f becomes greater than 9,999. 

Next an exponent ^ = (l/lO+J/2540. 62832) is formed. Next the num- 

r ^ r 

her 0 = y is formed. Next the number = (fractional part of o ) is 
obtained. Finally the relative error is determined as x r, where r is 
the preselected maximum relative error. Then 

e, . =* + 7) X r X a. . 
ij - i iJ 

The algebraic sign is given by the first digit in the decimal ex- 
pression of . If this is even, the sign is positive; if this is odd, 

the sign is negative. 

A special treatment is used if a. . = 0, for which e. . = + >) , the 
sign being determined as before. This has obvious disadvantages if the 
average value of the elements in A is small, and it is suggested that 
future work make an appropriate modification to the present procedure. 

In retrospect it appears that perhaps a more truly random variation 
could have been attained by using ~ (fractional part of 100 cT ), 
say. Furthermore, it would have been desirable to incorporate such ran- 
domizing influences as the date, or local time, or job sequence which 
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should become available internally to the computer, but which were not 
available at the time this work was done. 

For convenience in comparing values of C* with C, to see how good 
the recovered matrix C* was, most of the compliance matrices C used were 
diagonal matrices. However, there is no loss in generality as may be 
seen from the following argument. 

Starting with the matrix equation CU = V, with diagonal C, and pre- 
multiplying this equation by an orthogonal matrix M, it can be written as 

MCM*^ = MV (19) 

T 

where the insertion of M M in equation (19) between C and U does not al- 

T 

ter the equality, due to the fact that M M = I. 

Now let 

= C 
MU = U 
MV = V 

equation (19) can be written now as 

CU = V (20) 

Equation (20) has the same form as original matrix equation CU = V, 
but with the difference that starting with C (a diagonal matrix), the 
new symmetrical matrix C is no longer a diagonal matrix, and instead all 
its elements are generally non-zero. 

Thus a general problem CU = V could be written in the particular 
form CU = V, with diagonal C, if only the orthogonal matrix M, which dia- 
gonalizes C were known. C = M*^CM, U = M^^U, and V = M*^V. 

The particular validity of this argument was verified by two types 
of numerical experiment. In the first type, nondiagonal C*s were used, 
errors were introduced in U and V = CU, and the recovered matrix C* was 
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successfully compared with the original C* In the second type, an ortho- 
gonal matrix M was generated, and matrices C, U, and V were formed using 

M, diagonal C, U, and V. Then random errors were introduced in both sets 

- T 

of data (U, V and U,V) and C* and C* were obtained. Comparison of M C*M 

with C was quite satisfactory. 

From theory the minimum number of independent pairs of vectors F and 
D needed to solve the problem is 4. This is because in the compliance 
matrix C there are 21 unknowns, and 21 independent equations must be form- 
ed. Thus would require four pairs of 6-vectors F and D since the notion 
of 3% pairs is excluded. However, the particular methods employed in the 
programs seem to demand a larger number. The minimum number of pairs of 
vectors needed to obtain **good** results was six. 

In one case, a problem with 24 pairs of vectors was solved twenty 
times, using a maximum random error of ten percent. The errors were in- 
troduced in a different manner for each of the solutions. From the twenty 
compliance matrices obtained in this way, an average compliance matrix 
was obtained. This matrix looks ’Very good”, when compared element by 
element with the original compliance matrix. This indicates that if the 
value of n is so large as to overload computer memory or as to prove in- 
convenient for loading data, the U and V matrices may be conveniently 
partitioned, a number of ’’smaller problems” dealt with, and the answers 
averaged. One could regard the case described as forming a single pro- 
blem with n = 480 which was actually treated as the average of the re- 
sults of twenty problems, each having n = 24. 
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SECTION 8 



CONCLUSIONS AND RECOMMENDATIONS 

There does not seem to be any direct way of testing the validity 
of the theory and the computer programs which implement it other than by 
taking actual examples. These examples were constructed by assuming 
values for matrices C and U and generating matrix V = CU. Then, in most 
cases, errors were applied so as to get modified matrices U* and V^. 

Then, from these matrices, considered as input data to the program, the 
various procedures were used to recover an optimal value, C^, which is 
then compared with the original C. No useful or convenient method has 
been found for assessing the accuracy of the recovery other than by com- 
paring term by term, and, to simplify this task, in most cases C was 
assumed to be a diagonal matrix. It was shown in Section 7 that there is 
no loss in generality in doing this. 

One of the things that might be expected to influence the ease and 
accuracy of the recovery is the relation between the eigenvectors of the 
matrix C. Tests were conducted, at the one extreme, with diagonal C 
having identical diagonal elements, and at the other extreme, with dia- 
gonal matrices having diagonal elements differing by a factor of 3125. 

In all, fairly extensive numerical work was investigated and various 
maximum percentage errors were introduced. The output produced as a 
result of these studies is much too volumnious to be reported in detail 
and it has not been found possible to reduce the volume of information 
in a satisfactory quantitative way. Accordingly, it is necessary to 
restrict remarks to those of a general nature which will be made in the 
following. In general, it may be concluded that the theory and the pro- 
gram have been proved to be accurate and effective, and in combination 
they provide a useful tool for analysing the data from certain types of 
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structural and mechanical experiments. 

The following matrices have been used as compliance matrices C. 
from which the matrices O' v/ere recovered. 



( 1 ) 



(3) 



9. 


0 


0 


0 


0 


0 




1. 


0 


0 


0 


0 


0 


0 


00 


0 


0 


0 


0 




0 


2. 


0 


0 


0 


0 


0 


0 


5. 


0 


0 


0 


(2) 


0 


0 


3. 


0 


0 


0 


0 


0 


0 


12. 


0 


0 


0 


0 


0 


4. 


0 


0 


0 


0 


0 


0 


8. 


0 




0 


0 


0 


0 


5. 


0 




0 


0 


0 


0 


27. 




_0 


0 


0 


0 


0 


6. 



1. 


0 


0 


0 


0 


0~ 




~7. 


0 


0 


0 


0 


o“ 


0 


10. 


0 


0 


0 


0 




0 


7. 


0 


0 


0 


0 


0 


0 


100. 


0 


0 


0 


(4) 


0 


0 


7. 


0 


0 


0 


0 


0 


0 


100. 


0 


0 


0 


0 


0 


7. 


0 


0 


0 


0 


0 


0 


10. 


0 




0 


0 


0 


0 


7. 


0 


0 


0 


0 


0 


0 


1. 




0 


0 


0 


0 


0 


7. 



(5) 



5. 0 1. 0 

0 5. 0 0 

1. 0 5. 0 

0 0 0 5. 

0 2 . 0 0 
1 . 0 1 . 0 



0 

2 . 

0 

0 

5. 

0 




5. 




0 0 0 0 1 . 

2 . 0 0 1 . 0 

0 3. 1. 0 0 

0 1. 4. 1. 0 

1. 0 1. 5. 1. 

0 0 0 1 . 6 . 



(7) 



1 . 0 0 0 0 0 

0 5. 0 0 0 0 

0 0 25. 0 0 0 

0 0 0 125. 0 0 
0 0 0 0 625. 0 
00000 3125. 



OUTPUT No. 1 Gives the desired compliance matrix that satisfies 

criterion one. With this procedure the types of 
matrices numbered 1, 4, 5, and 6, have been recovered 
quite successfully using up to 35% maximum random 
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errors. Matrices numbered 3 and 7 were recovered with 



up to 57o maximum random errors, and matrix 2 was 
recovered with errors up to 25%. (Where no errors 
were introduced, each original was accurately recover- 
ed) . 

OUTPUT No. 2 Gives the desired compliance matrix that satisfies 

criterion two. All the matrices were recovered same 
as the originals when no errors were involved. Matri- 
ces numbered 1, 2, 3, 4, and 5 have been recovered 
quite successfully using up to 35% random errors, 
and matrices numbered 6 and 7 with up to 25% random 
errors . 

OUTPUT Nos. 3A,3B, Gives the desired compliance matrix that satisfies 
3C 

criterion three, using the procedure that obtains the 
square root of a non-s 3 nnmetric matrix. All the matri- 
ces were recovered same as the originals when no errors 
were involved. With matrices numbered 1 to 6 this 
procedure has been worked successfully with up to 35% 
random errors, and with matrix numbered 7 with up to 
20% random error. 

OUTPUT No. 3D Gives the desired compliance matrix that satisfies 

criterion three, using the procedure that takes square 
roots of two different symmetric matrices. Matrix 
numbered 1 is the only one that has not been recovered 
same as the original when no errors were involved. 

In the light of other results obtained, this failure 
to recover matrix 1 to a very high degree of accuracy 
is unexpected and it has not been satisfactorily 
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explained. Matrices numbered 1 to 6 using up to 35% 
random error, have been recovered quite successfully, 
and matrix numbered 7 with up to 257o. 

OUTPUT No. 4 Gives the desired compliance matrix that satisfies 

criterion four. Using any value for BETA, when no 
errors were involved, all the matrices were recover- 
ed same as the originals. When the error was increased, 
better results were obtained for small BETA*s, The 
same remarks given for OUTPUT No. 1, applies to this 
output. 

To have a better idea how these results look, see Appendix C, where 
the computer solutions for a particular problem with zero and ten percent 
random errors, are given. 

The work presented here is just the beginning and not the last word 
in this new area; more work and investigation should be done. Some sug- 
gestions are given below: 

1) Solution of equations that satisfy criteria 1, 2, and 4, 
can be investigated using the procedures given by Bickley and McNamee 
[4], or possible modifications to take advantage of the fact that matrices 
A and B are equal. 

2) Investigate how large maximum random error percentages one 
can give to matrices U and V. Give to these matrices different maximum 
percentages of errors, e.g., to matrix U, 5%, and to matrix V, 15%. 

3) There exists a fifth criterion given by Brock [1], that has 
not been investigated yet. 

4) Devise a new iterative procedure for the equation that 
satisfies criterion four. 
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5) Devise a convenient criterion which measures the ’’overall** 



accuracy of recovery of a matrix C* At present, one merely compares 
element by element. 

6) Employ advanced and sophisticated statistical theory for the 
analysis of errors and the accuracy of recovery. This would seem to 
imply constructing an adequate set of definition and the use of probabil- 
ity spaces, etc. The present analysis, in which general conclusions are 
attempted from specific cases, leaves much to be desired. 

As things stand now, the engineering analyst has, in this theory 
and program implementation, a four edged tool which can be used for the 
analysis of experimental data. The decision of which criterion yields 
the ’’best” result is one which must be made on the basis of experience 
with the physical problem. 
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APPENDIX A 



INSTRUCTIONS FOR USE AND LISTING OF PROGRAMS 
A-1 Preparation of Data Deck and Final Assembly of Programs. 

The preparation of the data deck is the same for programs CECI 1 
and CECI 2, but has some alterations depending on whether MODE = 1 or 
MODE = 2. 



In what follows the instructions of how to prepare the data deck is 
given for both types of modes. 

The first card is bbbbbbbbNNbbblbb where b = blank, the 6 is 

the number of rows, NN is the number of columns, and 1 means first matrix 
to be read. When MODE = 1, NN = 06, is the number of columns in matrix 

C; when MODE = 2, NN varies from 06 to 50, and is the number of columns 

in matrix U. (Thus U is limited to a maximum of 50 columns). 

Following this card, comes the elements of matrix C or U as the 
case may be, punched in five entries to a card, in FORMAT (5F14.6). 
Entries for the first row are thus entered until all elements for the row 
have been punched. Then, starting with a new card, entries for the sec- 
ond row are punched, and so on, until the entries for the sixth row have 
been entered. Each row starts with a new card. 

Card bbbbbbb6NNbbb2bb , follows. Then the entries for the 

matrix U (if MODE = 1) or matrix V (if MODE = 2) are punched in the same 

manner as before. 

Example: Punch cards as in Figure A-1, for the following 6x6 diag- 

onal matrix as C matrix. 



C 



1 . 

0 

0 

0 

0 

0 



0 

2 . 

0 

0 

0 

0 



0 

0 

3. 

0 

0 

0 



0 

0 

0 

4. 

0 

0 



0 

0 

0 

0 

5. 

0 



0 

0 

0 

0 

0 

6 . 



44 




Note that element C(l,l) is punched between fields 1 to 14, element 
C(2,2) between fields 15 to 28, and so on. Note too that each row has 
two cards because there are six elements for a row, and in each card 
only five consecutive elements can be punched. 

If more than one set of data is used, the second set comes after 
the last card of the first set, starting again with the first card 

bbbbbbbNNbbblbb 

and so on . 

In order to stop the program, the final card of the data deck is 

punched as follows bbbbbbb0b6bbb6bb 

The complete program deck includes; control cards, main program, 
subroutine package, and data deck. The assembly of these components is 
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shown in Figure A-2, where the first control card is from the CECI 1 
program. To use program CECI 2, change the word CECI 1 to CECI 2. 




FIGURE A-2 
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A-2 Listing of Program CECI 1 



The cards of program CECI 1 are numbered from 0000 to 0159. 



o ivj nJ- iT "0 r- oc O' o rvi ro if\ vO OD Q <vj lo ^ 0 ^ O' o (M rr 1 St iTi o ® O' o fvi m IT <o r- 

o o o O o o O O o o ^ 1 -^ ^ ^ ^ M <N rvj c\i rvi <\i cvi <vj c\i fvi fn ro fTi rn rn rn fn fo ^ >t 

oooooooooooooooooooooooooooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 



•» w 
^0< 



irsfsj 

wfvj 

QC^ 

h-«J 

•►o 

-"O 

o ► 

w 
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A-5 Listing of Progreun CECI 2. 

The cards of program CECI 2 are numbered from 00000 to 01570. 
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APPENDIX B 



SUBROUTINE PROGRAMS 
B-1 Description of Subroutines. 

SUBROUTINE ADDEM (A,B,C,D,M,L) provides the result of A = B_fDC, 
where D is a scalar constant, and A, B, and C are double precision matri- 
ces dimensioned (M,L), where M and L are not less than 3. With D = 1.0, 

A « B-tC; with D = -1.0, A = B - C. 

The matrices B and C from the calling program must have identical 
order, otherwise a message is printed and the execution is halted. 

SUBROUTINE AINVM (A,B,C,D,L,N) forms the inverse A = B"' of input 
matrix B. Matrices A, B, C, and D are double precision, with A and B 
dimensioned (L,N). Matrices C and D, which are used as scratch matrices, 
have dimension (N,N), where L = N+1, and N is not less than 3. 

A message is printed if the input matrix B is non-square and the 
execution is halted. 

SUBROUTINE BIGCM (C0L,Q) forms the 21-element vector COL from the 
6x6 matric Q, according to the scheme described in Section 3; if the 
input matrix Q is non- square a message is printed and the execution is 
halted. 

SUBROUTINE BIGSM (BIG,P) forms the 21 x 21 matrix BIG, from the 6x6 
matrix P according to the scheme described in Section 3; if the input 
matrix P is non-square a message is printed and the execution is halted. 

SUBROUTINE BLANCA (CO,BE,N) forms matrix CO from matrix BE by intro- 
ducing random errors in the latter according to the method described in 
Section 7. The arrays CO and BE are dimensioned (7,N). The maximum 
relative error, denoted by the symbol r in the exposition in Section 7, 
is governed by the computer variable PERC, which is one tenth the relative 
error. If PERC = l.D-2 = .0100..., then r = .100..., and a maximum of 
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10% error is introduced. PERC may be changed by changing card number 
0313 in the subroutine. 

SUBROUTINE JACVAT (B,N,N0YES, RR,SS,EIRU,NDIM) gives the eigenvalues 

and eigenvectors of the symmetric matrix B. If B is non-symmetrical the 

T 

subroutine first makes B symmetric by forming %B+^B . RR is the spectral 
matrix that contains the eigenvalues of B. SS is the modal matrix of B 
that contains the eigenvectors. EIRU is the square root of the elements 
in the diagonal matrix RR. N is the dimension of B. NOYES is 0 or 1; 
if 0, no eigenvalues will be recovered; if 1, the eigenvalues will be 
recovered in matrix RR. In the application made in this thesis, NDIM is 
the same as N. All the matrices and variables are in double precision. 

SUBROUTINE MAUSS3 (N,EP,A,X,KER) is a double precision FORTRAN IV 
language subroutine which inverts a symmetric matrix. A is the matrix 
to be inverted, X is its inverse. Both have dimension of (N,N) . EP is 
a tolerable error, taken as l.D-20 in the present work. KER is an output 
flag having the value 1, if all has gone well, or 2, if there has been 
any difficulty of execution. 

SUBROUTINE OUTPU (R,B,U,V,T1,T2,T3,T5,LL,L) has inputs R, B, U, V, 
where R is the desired compliance matrix . B is the desired stiffness 
matrix, and U, V are the data for the problem (such that RU = V, approxi- 
mately) . It prints the matrix R, and computes and prints the three 
criterion numbers o(^ (i = 1,2,3). Matrices T1-T5 that appear in the 
calling statement are scratch matrices. LL and L are dummy variables 
used in the dimension of the matrices and are related to the number n of 
vector pairs used. 

SUBROUTINE PRINM (A,L) causes the (6,L) matrix A to be printed, 
where L is not less than 3. The first line printed contains the fol- 
lowing information: MATRIX NUMBER XX, YY ROWS, ZZ COLUMNS. This is 
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followed by a print of the matrix, line by line, in FORMAT (IX, 8D12.4). 

SUBROUTINE RECIP (A,B) A and B are diagonal matrices (6,6) and A = 

B where B is the input matrix. 

SUBROUTINE REDEM (A,IDENT) reads into the computer memory from the 
data deck which is punched as described in Appendix A. A is the "name” 
of the matrix and IDENT is an identifying integer to assure that the 
card-group which is read is indeed the one desired. 

SUBROUTINE SQRERT (A,B,C) solves the matrix equation that satisfies 

criterion 3, using the procedure that twice takes the square root of a 

T T 

s)mimetric matrix. B UU , C “ W are the inputs, and A is the desired 
compliance matrix. 

SUBROUTINE SQROOT (A,CC,G) takes the square root of a non-s 3 Tnmetric 
matrix by iteration procedure. Matrices CC and G are inputs. CC is the 
approximate square root; G is the matrix whose square root is desired, 
and A is the desired square root. 

SUBROUTINE SQUAM (A,B) forms a 6 by 6 matrix A from the 21-element 
column vector B, arranging the results properly as implied in Section 3. 

SUBROUTINE TRACM (T,A) foms the trace T of the 6 by 6 matrix A. 

SUBROUTINE TRANM (A,B,LL,L) forms the transpose of matrix B(7,L), 
obtaining matrix A(LL,6), where L is no less than 3, and LL = L + 1. 

SUBROUTINE ULTMA (A,B,C,MM,M,L,LL) forms the matrix product A = BC. 
A(MM,M) = B(MM,L) x C(LL,M), where MM, M,L,LL, can be any numbers not less 
than 3 for which the matrices are conformable for multiplication. A 
message is printed if the matrices B and C are not conformable and execu- 
tion is halted. 

SUBROUTINE ULTMB (A,B,C,LL,L) multiplies the square matrix B by a 
column matrix C, to obtain the column matrix A. A(LL) = B(LL,L) x C(LL), 
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where L is any number not less than 3, and LL = L + 1. 

SUBROUTINE OUTPUT (R,B,U,V,T1,T2,T3,T5,LL,L,BETA) is the same as 
SUBROUTINE OUTPU, The only differences are the inclusion of parameter 
BETA in the calling statement, and the printing of four criterion numbers 
instead of three, (i = 1,2, 3, 4). 

B-2 Subroutine Listing 

CARDS NUMBERED 



SUBROUTINE 


FROM 


TO 


ADDEM 


0160 


0182 


AINVM 


0183 


0208 


BIGCM 


0209 


0233 


BIGSM 


0234 


0288 


BLANCA 


0289 


0336 


JACVAT 


0337 


0481 


MAUSS3 


0482 


0532 


OUTPU 


0533 


0559 


PRINM 


0560 


0576 


RBCIP 


0577 


0594 


REDEM 


0595 


0616 


SQRERT 


0617 


0638 


SQROOT 


0639 


0692 


SQUAM 


0693 


0727 


TRACM 


0728 


0744 


TRANM 


0745 


0763 


ULTMA 


0764 


0791 


ULTMB 


0792 


0814 


OUTPUT 


20000 


20360 
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APPENDIX C. 



SAMPLE PROBLEM SOLVED BY COMPUTER 

C-1 Sample Problem Using Program CECI 1 with 07© and 107® Maximum Random 
Errors. 



The nomenclature used in the output of the program is as follows: 



MODE = 1 
Matrix P 
Matrix BE 
Matrix EB 
Matrix U 
Matrix V 
Output No. 1 

Output No. 2 

Output No. 3A 

Output No. 3B 

Output No. 3C 
Output No. 3D 



Data will be generated in the program. 

Original matrix C as read from cards. 

Matrix U as read from cards. 

Matrix V from V = CU. 

Matrix BE with random errors introduced. 

Matrix EB with random errors introduced. 

The desired compliance matrix which satisfies 
criterion one. 

The desired compliance matrix which satisfies 
criterion two. 

The desired compliance matrix which satisfies 
criterion three, using the iterative procedure for 
taking the square root of the non-symmetric matrix 
wW. 

The desired compliance matrix which satisfies 
criterion three, using the iterative procedure of 
taking the square root of the non-symmetric matrix 

uu^w*. 

The desired compliance matrix which satisfies 
criterion three. It is the average of outputs 
3 A and 3B. 

The desired compliance matrix which satisfies 
criterion three, using the procedure that twice 
takes square roots of symmetric matrices. 



The printed output for the problem with zero error appears on pages 



77 to 82 inclusive; that for the problem with 10% maximum error ap- 



pears on pages 83 to 88 inclusive. 
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C-2 Sample Problem Usint? Program CECI 2 with 0% and 10% Maximum Random 
Errors . 

The nomenclature used in the output of the program is as follows: 



MODEsl 


Data will be generated in the program 


BETA 0.8 


(The sequence used is .1,.2,...,.8.) 


Matrix P 


Original matrix C as read from cards. 


Matrix BE 


Matrix U as read from cards. 


Matrix HIl 


Matrix V fraa V ~ CU. 


Matrix U 


Matrix BE with random errors introduced. 


Matrix V 


Matrix EB with random errors introduced. 



Output No. 4 The desired corapliance matrix which satisfies 
criterion four. 

Tne printed output for the problem with zero error appears on 
pages 90 to 93 inclusive; that for the problem with 10% maximum 
error appears on pages 94 to 97 inclusive. 
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MATRIX EB FCLLCUS 



ro rn rn fT, CO st fT' ro (NJ 

ooooooocoooooooooo 

ooooocoooooccoccoo 
00 %c o 00 o 00 o o in o o CO o o '■o o in 
^.-^lnf^o^-lno^'0(^JO^oaoc^^'0 
m sO O rein fn >o m ro h* ^ >0 o — < nj nj 
^ in ^ O' ^ m O' re vO nj f-< in in 



oooooooooooooooooo 
II I III 

-^rernrerofe rvj cnj re corn cere rn^rn < 4 * 
oooooo oooocoooooo 



I fsj fv nj ^ f\j r\i f\j f\j f\' fsj fsj fvj »-♦ r\j r\j tNj r\i 
OCOOOO OOOOOO oooooo 

OOQOOOOOOOOOOOOOOQ 

oooooooooooooooooo 
esj o ^ o ^ o o m o o c^ o o 00 o o m 
O <M in nj in eg e* rvj >0 00 in ce m ' 4' 00 rvj 
in ^ in eg CO 00 re %c re e- 00 eg ^ rvj 



oooooooooooooooooo 

II I III 

oegegegegeg eg ^ eg eg nj eg eg rg eg eg eg 
oooooo O OO OOOOOO oo 



oooooo ooooooooooo 
oore^g-oo oooovoooooreo-^ 
oego>^^>n-^ oooeega>j-oo^orere 
reinme-oeg insO>nrnc 30 0 'CO’ 4 'ooo 
>0 eg re re re %o o h- ^ e- e- in eg h* 



oooooo ooooooooooo 
oooooo ooooooooooo 
ooe-oooin oooooooo-^ooeg 
oooNCoe-’^’ inego^—unre'-^rea'a^co 
e* rg re eg re O re re nO > 0 eg sO oc eg re 



oooooooooooooooooo 
II III 



oooooooooooooooooo 
II III 



re re re re re re eg re re re rervirere"^ 

o ooooooooooo ooooo 



eg eg eg eg eg eg eg eg eg eg eg eg*-^egegeg 

ooooooooooo o ooooo 



Q ooooooooooo ooooo 

h-Oino>OegoOOOOO OOOO^ 



oooo-^O'Oininooegego ocoinineg 
— -O^Qcegref^Of<*)^re>o egcx)4'00 
eg ^ nj ^ re eg sO ^ O — » >0 O' 



oooooooooooo ooooo 
oooo oooooo oo ooooo 
reoinino-^oooooo oooooo 
Ntrgreoh-ininoOi-»>oin ing^inine- 
eg*-^egi-i>oe-re%o>t‘in re re 1-4 re 



oooooooooooooooooo 
I I I I I I 

re rere ^ ^ «-irere re re reeg egeg ^ ^ 
oooooooooooooooooo 



oooooooooooooooooo 
I I I I I I 

eg eg ej eg eg o eg eg eg eg eg eg eg f-i eg eg eg 
OOOOOOOOOOOOOOOOOO 



OOQOOOOOOOOOOCOOOO 
oO'4'Noegoooinoooooo*-^a'>o 
in ^ > 4 - eg o o O O rg sO ® O o O vO re 
o e- re r* O ^ eg O' o e- in o ® eg in eg in 
00 f-i C" O ^ eg eg >}• ® ^4- ^ ^ »-4 



OOOOOOOOOOOOOOOOOO 
oooooooooooooooooo 
O 0 ^ 000000^0 OmOOO'O OCT 

O*' in C^ oeg ® in re ^ rein nO 

h- ^ O' e- in %0 rsj O' ® eg eg in in in in 



oooooooooooooooooo 
I II II I 



oooooooooooooooooo 
I II II I 



re re eg > 4 - re reeg re eg rerererere' 4 *re% 4 ' 
000000000 00000000 



(M eg eg eg eg eg eg eg eg eg eg eg eg eg —4 eg 
000000000 00000000 



00 

Z 

z 

o 

o 

o 

' 4 - 

eg 

•k 

I/) 

JC 

o 

QC 

-O 



OQO OOOOOO 00000000 
000 >co 0000 OOOOO® oeg 
eg O' >0 ' 4 - O' O mo eg>ogD®®egsOm 
'4'm®mrecrmgDm mo'mosoeg-^o 
re s4* re ^ eg O' O' eg e- o eg >o eg nJ- eg eg eg 



oooooooooooooooooo 
III I 

re re re >4* '4' eg re re re re re re re re re >4- 
o 0000000000000000 

a 000 QO oooooo aaoQa 
O ^ 00 ceo OOOOOO 00 ego® 
>ro-^r--oomoeg>reg^oomeg 
h- > 4 '®r^-^e*mm®a'omO' 4 -*-^rer^ 
'-yO'Orei— <•— <«-ym'*^'-g>4'Com»^'4'NO**^'^ 



z 

X 

o 

o 

o 

> 4 ' 

eg 

<>o 

jg. 

O 

esC 

>0 



000000000 00000000 
000000000 00000000 

ooooomooo oooommoo 
Op-^ree-remerrem ^®eg-^remo® 
rem>4'®»-^m*-^m^oegmremre^4' ®n- 



oooooooooooooooooo 
III I 

eg eg eg eg eg eg rg eg eg eg eg eg eg eg eg -4 eg 
o 0000000000000000 

Q QQOOQQGQQQOOOOOQ 
o 0000000000000000 
o ®oomooo oooooo® 00 
O' ^omm>4'*-4— »ie\*-^r*o'remego'4' 
O n- eg O' ® re re ^ > 4 - ® ® m eg m ® 



oooooooooooooooooo 
I I I I I III 

•• 00 

re re re re re re re re re re re re eg re re *-4 re re re • • 

oooooooooooooooooo o o 



oooooooooooooooooo 
I I I I I III 

•k 

' 4 ' eg rvj eg eg eg eg eg eg eg eg eg »- 4 rsj ego eg eg eg 

oooooooooooooooooo 



QC 

uu 

Oi 

X 

3 

z 

X 



QC 

»— 

<I 

X 



OOOOOOQOQOQOOOOOOQ 
0®'4*0egeg000000000000 
® m >4* o ® m m o ® eg o ® eg o >4* e- re 
ego'®-^®> 4 'rer^p-^ 00 '> 4 '> 4 '*-^>t®®® 
cx> h- ^ ® ^ ® ^eg ^4^ >4* eg >4* re ® ® m h- 



oooooooooooooooooo 
I II I I I I I 

re re re ^ re re re re re re re re eg re ^ 4 - > 4 * 
oooooooooooooooooo 

000000 a oooooo ooooo 
oooor^® 00000000® o^ re 
re®oegn-^omm®®egooreegoeg 
N»®m®r^®— 4 m> 4 '®®o^o®mm®re 
® eg >4* ^ O' eg rg re e- — 4 eg ^ re 



OOOOOOOOOOOOOOOOOO 

II III I I 



II 


M 


to 








:s 




< 


<l 


0 


QC 


0 


0 


-j 


Ui 


z 


Z 




OJ 


<[ 


< 


0 


X 




a 


UL 


3 


cu 


CO 


3 


Z 


z 


z 




X 






X 


►-4 






»-4 


QC 


QC 


QC 


QC 


K 


0 


0 


H* 


< 


QC 


QC 


< 


X 


QC 


QC 


X 




Ui 


Ui 







OOOOOOOOOOOOOOaOOQ 

oooooooooooooooooo 
O ' 4 * ®00' O' 000 00 00 000 00 
eg ® ^ m ® m e- m eg s 4 - o ® O O eg O' 

O' ® m k4* eg re eg re >4* re >4' rg m re ® re eg eg 



oooooooooooooooooo 
I II I I I I I 

eg rvj o ^ eg eg eg eg eg eg eg rg eg eg eg eg eg 
OOOOOOOOOOOOOOOOOO 

ooaaaaoaaoaaooooQa 
oooooooooooooooooo 
oooof^r^ooooooooegreoo 
e- eg o o ® >r nj O' ^ > 4 ^ ® m o >r ® h* O' 
O'remoa'm%4*m®^®'-^egeg4'm®-^ 



oooooooooooooooooo 

II III I I 



91 



SNiwmDD 9 ‘9 yagwRM 



ooccoooooooooooooo 

OQQOCOQOQOOOOQOCOC 
00 O CO C C3D o C in o O cc O C Nt >c o IT. 

^ ^ iTy o m o ^ 'O fNj o ^ o cc st o 

m vO o in ro so rn ro h - \0 o osi O' CO 



ocooocoooooooooooo 

II I III 

•—•romrcrcfc (\ic\jrcrOfCrorOfc>^'(<'. 

ocoooo oooocoooooo 



QOCCCCi oooocccccco 
ooc^^oo ocoo'ccocomor-^ 
Of\JO>tvO*-< OOOOCOO'N^'OONt orr'rc 
rrinroc-oco m'OsOrccoO'CO^cx5 O 
NOcomrcfC^oor^^^c-c-— ^^ir.cviC- ^ 



oocoooooooocoooooc 

II I I II 

rcrorcrc>^'>^rcrcrsjrcrcrc rocsircrcvt 
OO COO O OCOOOO oooco 



OCCOOOOQQOCQ OQQQO 
C'OmOsCfNJOOOOOO ooco^ 



QOoorHO^oininoococjo ocoimncj 
r-^o»-'xcsifor^oro-^rcsO cnjcd'^'OO 
f\l f\l ^ rc >0 ''O O' 'O O* 



cooooooooooooooooo 

I i I I I I 

rOf<'.rnNtNt^fCrnfn'^rCfncsi(NJoo>cf>t >t 

oocoooooooooooocoo 

OOCCOOOOQQOCOOQOQO 

OC>3- orvjcoomoooooo^o nO 
^ LT t >t Csi O O O O CO sC CO O o o vO C". 

O' r^mr^ >^*0'— ‘fMOor^mocofNimcsj m 

CO sC ^ O rC r— < »—t f\J f\| CO < *— • 



OOOOOOOOOOOOOOOOOO 

I II II I 



r<*iror\}>trCrcf\ir<*NCsi 

ooooooooo 



c^rCrcror^v4'fn>^' 

oooooooo 

OCOCOCQC 

OOOCOCOOCvj 

cj>cocx:co(\isom 
incrincjsOfNjr-^o 
f\j f\j ^4- f\i ^ r\j f\j 



ccccocooo 

OOOvOOOOOO 
cja c-vONf ooino 
NTincOLoroo^LnsOm 
(<'. vj- rc f\j cc f\j C* O 



2- 



OOOOOOOOOOOOOOOOOO 

III t 



rn ^Rr^v^>j-rco^cornrnrC'^rnrn(<'.rn <r q 

CJ oooooooooooooooo 



o cjooaooaaooocjcooo 

0 >r cj-occooooooooonjox) 

— < xrcj--^f^oomof^J'rf^J^TOOl^v(^J 
0- Nro^-^r-u“Mnxc^Olr^o^^^mr- 

«Ln^r-i>^com*-^-4- vO^ ^ o 

cr cr 

OOOOOOOOOOOOOCJC^OOO • • 

1 I I I I III C' o 

rcrcr<^rnrn(^.rnc^c^c^rcrvjrcr<^r-^'^rcro 
c:5 0oocj OOOOOOOO ooooo 



OOaaQOOOQQOOQOQQQQ 

o -O O CsJ f\j o O o o o o o o o o o o 

couC'^ONr'Omir>ocof\joco(\jo>rr-m 

f\icrso^x>Nfroc«*HOcr'>tsr'-<N3’NO<Doo 

coc-srxN^NO-^^Csi>j*s^'r\jNrroNoaomr- 

• • • • II II 

OOOOOOOOOOOOOOOOOO 

I II I I I I I LU LU 

r X 

rc rc on (n rc rc rn ro ro rn rn f\j nJ- nT vi* ^ ►-• 

OOOOOOOOOOOOOOOOOO K- h- 

OOQQOOOOOQQOQOQOao O O 

ooooc-^oooooooosoocrro uj uj 

^xOf\ir- '4'OiruoxxfMOOrorsjocg c<o ui 

^xir\>cc-uD»-'j^'>rNC>ocro>CLn«rxc'j a_ cl 

XCJ>f^'--<OCJCsjrOr^f^^f\|r^fr.r-i^r-^ < <J 



OOOOOOOOOOOOOOOOOO UJ LU 

II III I I 



92 



NUMBER OF ITERACTIONS = C.80CC0CCC 



06SIRFD MATRIX IS 



^ C. 

1 i I I I 

oocooo 
lr^oou^f^J^to 
^Q0Q0 00^tO 

>0000 Qoor^ 

• ••••• 

oooooo 

I i 

iTkinin ^ —I m 

Mil I 
ocoooc 
^00 

00 ro ITI O" C ^ 
oostocroo 
cr, 0> ro ^ 00 » 



CO 

q: 

o 

a: 

cc 

OJ 

uu 

Qc 

< 

o 

o 

CO 



LU 

cn 



> 

CO 

LU 

X 



oooooo 
I I 

NCst 

^ 

lit It 



KT. 


COOOOQ 






in 


'J' 


rv! 




2" 


sOCT ooof\j 




CNJ 


(NJ 


nj 


(N 




X 


f\j>r vooo'oo 




1 


1 


I 


1 




O 


mroo(\(0^oo 




a 


Q 


O 


O 




-J 


OL' »— < 1-^ 00 




o 


O' 


n 






o 






GO 




r- 






o 


oooooo 






>r 


00 


o 






1 1 






un 




o 










O' 


in 


ro 


o 






m nJ* •-4 tn 






in 


ou 


o 




«» 


»— *< <— 4 C , < »— H »— H 




« 


• 


• 


• 




W) 


It III 




r-^ 


ro 










CJQOOQO 














U 


u'CT'Otrrom 














0^ 


cn CO rD nO iTk 00 
















r^ooox 




C^ 


CO 


m 


cn 




o 


r^rncn^ro>4* 




-J 








O" 








< 


< 


< 


c 


O' 




OOOOOO 




O 


o 


ZD 


o 


• 




I 1 




o 


Of 


Of 


o 


o 


*• 






aj 


LU 


LJ 


LU 






LT. (\j nT nT O 




















p-4 


Osj 


f<) 


sT 






1 1 1 1 I 

OQOOOQ 
















coao'sor^ao 


O 














cooo ^r (^°o 


o 














moor^f<>>r o 




a: 


QC 


QC 


QC 






r—4 r—I fO r—H *>*^ 


o 


aj 


LL’ 


aj 


LU 




X 


•«•••• 


o 


cC 


CD 


CD 


CD 


II 


LU 


OOOOOO 


o 


X 


X 


X 


X 




OJ 


1 ( 1 


o 


o 


o 


o 


o 


LU 


X 




o 




x: 


x: 


X 


X 


o 


•— ‘incn'CLrv'«r 


o 










— « 






00 






X 


X 






1 I 1 1 1 


• 


u 


o 


u 


u 




X 


oaoaoQ 


o 




t— 4 




►-< 


o 




c)ooir\r*-vom 


li 


QC 


QC 


C£ 


QC 


u 


00 


ococnr\jc30>o 




LU 


UJ 


LU 


UJ 


c^ 


h- 


orn^^'^f^-ooso 


<. 


h- 


h- 


►- 


I— 


a« 


< 


0> ^ CO rfi (Ni 


H- 


H- • 




*— < 




<t 


X 




OJ 


X 


QC 


oc 


QC 


-1 




oooooo 


CD 


o 


O 


o 


O 


UJ 



NO 

uc 



I I I I 1 



93 



TOTAL ELAPSEC TIME 



MATRIX P FOLLOVjS 



<N <\J <\J fV fVJ <N (NJ fV CM (\J ^ CsJ (\J fVJ (V 

ococcoooccc coocooo 



cccooocooooooooooo 
O O C C C C-> C C O C C C O O C C c. 
CM ^ o o o o fc> o o cr o o 00 o o 
O CM ITN fvj IT <N h- CM nO 00 «-• ITN cn m CD CM C^ 
IT IT CM CO 00 res Oh- C^^ h-ODCM^-^^CNJ--* 



COOOOOOOOOC'OOOOOOO 

II I III 

OCMfM(N<V'(M rvJr-4rMrMrM(\JCvt(\J<N(MCM 

oooc oo oooooooooce 



ocococ. 

ooocoo 

ooh-cooir. 

ocooof^'«t 

r^fMrorvj-Hmo 



CICCOCCOOCOO 
ooocooooooo 
c o cr o oc o o .-H o o (\ 
infMCT'-'incn^rc'.O^a'oo 
*~trofOvOO*-‘<\jvOoocMfn 



oooc oooooooooooooo 

II III 



CM 

c 



CM'fM<N<\JCM'CM(\jrM— »CM<NCM <\J *--• CM CM CM 

OCOOOOCOOOOO ocooc 



c 

o 

o 



OC COOCM 



OOOQOOOOOCCC 

oooococoocoo 
m o ir> tn o o o o cj o o 
st<\ 3 rooh-trLOOvC«-»'^ir. 
(\j^(\j,-i>^ 5 r*-cnvO'tir(^LO 



ooooo 
ooooo 
ooooco 
tTsOirur^h- 
O— *oorOi— ifTi 



oc ooc o 



oocooooocooceooooo 

I I I I I I 



— • CM CM <N CM CM O CM rj CM CM CM CM---' CM .-H CM CM oo 

c oocooooocooc ooococp 



c: 

c 

C- 

o 

COC-CCCC 



OOQCOOOOOCCCCCOOOC 
o ocooc ooc ooc ooc coo 
o o voo c o oo--*o o LT oo o nO c cr 

0 'tri' 4 'h-OOCM' 4 ’crLnr<^. — 'O'— 

Oh-'-HChr-LT vD<M croo CM (Min.-«LnLn>t IT 



ocooc.' o 



oocoooc^coc oeooocoo 

I II II I 



CM 

c 



<M(M*-'rMCv.:rMfMrvjf\j fMrvjCMCMtMrM^rM 

ooooooooo ocoooc:)oo 



o 


c 


oo 


I?' 


o 


z: 


2l 


o 


2. 


O 


CM 


O 


o 


OOC— <c c 


-J 


c 




o 


o 


OOOOOC 


o 


sD 










CM 




o 










JC 


c 




a 


c 


u 


O' 


o 


a: 




c 




sC 


ooinc.oc 


nO 



ecjooc'o 









•> 


—4 


CM 




CM 




C' 








Q 








o 








O 


on 






00 


3 






o o o o o 


o 




a 




o 


a: 


LU 


oc oooo 


o 


UJ 


OJ 




o 


OJ 


5: 




LI. 


2: 


o 


— • 




o 




o 


LU 


z 






cu 




X 


o 




X 




o 


X 


HM 




c 


»-4 


oc 


1— 


CL; 


a: 




< 


ooocoo 


k- 


< 


2* 




< 


3: 




oooooo 


2 





COOOQO.QOC 

ooooooooo 

cooooinoco 

oo^mr^romor^.io 



ooccoooo 

oc^oooooo 

OOOOlTMnOO 

0 <MmmLoro>rooh“ 



o ooooo ocooooooooco 

III I 



(M f\,(MCvjrMCvjCMrMrM^(M'f\i(MrMrvj^CM 

o ooooooooc c>oooooo 

O OQQOQO.QOCCOOQCQO 

o oooooooooooooooo 
o <?ooLT»cjoc:^oor oooncoo 

—lOh-cMOsomf^K^^vtNCvO^^ir.fMmsO 



oo oooooooooooooooo 

I I I I I 111 

CM CM tM (M (M <M CM fvj CM tM CM ^ CM O CM CM CM 

oocoocoooooooooooo 

QOQQQQOCQCiOOOaOQQO 

OO oooooooooooooooo 
O <3 oo (7^ oo O O O o O O O oo O 
CM 00 *— 'LO vO LO h- LO CM' ^ O O O O CM ^ O 
OocLONfrM''OCMmsrrONrfMunrooor<^CMCM 



oooooooooooooooooo 

I II I I I I I 

CM CM O —I CM CM CM CM <M CM CM CM CM -H CM CMJ CM <M 

oooooooooooooooooo 

OQOOQQQOOQQOQQQQOQ 
OOOOOOOOOOOOOOOOOO 
OOOOh-h-OOOOCOOOCMroOO 
h-(MOOQO> 4 'CM^cr'^>tsjir>o>r-oh-o 
O'coLOCT'OLn'j^LnNO'— •vO*~«cMCM>i'LnvO'^ 



OOOOOOOOOOOOOOOOOO 

II III I I 



94 



MATRIX EB FCt_LCV.,S 



OCOCCOOOOOCOOCOOOO 

000000000000000000 
cc>cooocooooircc^oooo>t'«ccir> 
-^ir O IT O ^ vO fNJ O C GO O' ^ 
ir> nO o o iT\ rr ' no m m r- —I %c o »-* o osj 

ococooooocooooooco 

II I ill 



ooooc e 



rv c\i rn ro ro rn O '4' 

cooc c ocoooo 



oocooo cocoooococo 
oor^>too oooovoooocn-io»-^ 
Cf\jc> 4 ^sOf-‘ ooocfvo^oD> 4 -orno 
ro LT fo o rvv* ir'NO'^mocONOO' 4 'OOO 
<) c\j rn rn rn sC o r- rM xf r- h- in f\j r- ^ 

OOOOOC'OOOCOOOOOCOO 
i i III 

rr<<\focn>t>4"f^fnrsjrnfr rn r^(vrri(n>4" 
C COOC OOC COCO oooco 

ooooooooocoo ocooo 
r^oinoocNJC ooooc oooo*-' 

ac<Tf-^(T'OirMrvooc\jfviO ooDirunrsj 

CNJCO^tOC 

c\J <•— i»^r<^rvj'0>4'v00 ^ \0 O*"^ 

cooooooc oooooooooo 
I I I I I I 

COOCOOC OOO OCOOO ooo 



OOOOOOOOOOQOOCOOOO 
C O > 4 - sOfV C O C in OC C OO 0-^0 vO 
— ‘in st (NJ OO O O c\j 00 O C st O O ro 

O' rr ^ 4 - C — » rvj O' c in o 00 f\j in c\j in 
oc >n pi — 4 *— < O ' ro > 4 " c\j nj * 4 “ oc > 4 " 

OOOCC C'OC oc OOOOOOOC 

I II II I 



CO c\j fVi fs . f-H fv c \ (\j <v (V fv r fM f -< <\j fv CM fv 

CCCCOCC C c OOC OC OCOC 

ooooooocooooooocoo 
cn (N. o a 00 4 -^ >t vO %t o r*- ^ 

c\j c> ^ in in o >n <v ^ sc in (\j IT in p^ 
fO on ^ CM sO p^ > 4 - nj >C ^ o >t nj o 

inm>tincMr^cm>cp^roP^QccMi-^> 4 -cMrvj 

OOOOOOOC C OOC OCOOOO 

II I lit 

O CM CM CM CM CM O (Ni CM PM rsj CM CM (M J CM' (M 
0000 C) 00000000000 C :0 

COOOOCCOOOCCCOCCCCi 
O' r-oir'«-^oc>oor)OG^cnNOooinsO«-^o 
^ in sO 00 nO 00 >C O (M P ^ ^ Nt CM oc 

ocroccMsorocmcrovO^cMCMp^ina'C 

P*-(MrOCM«~'(<',CMi-HCMN 4 'NC)r^^CMvCOOfM<n 

OCOC OOC CCOCOOC OOOO 

II III 

(M C\J (M CM f\‘ CM CM CM ^ CM CM (M o CM CM C\: CM 

O OCO OCOC COOOOOOC c o 
ooccococooooooocoo 

P -4 ^ vC O CMP^ IT o sO ir vC rr o O rn in CM 
rn PM cr >t fT in CM h- »-« sn >n s4- sT CM o P^ -H 

in pm st o ^ 00 — * o ^ CM Nt o o cn m cn in o 

fM^cMi— •p^r^rOvO> 4 ’infON 4 '> 4 "<— *orn^> 4 ' 

oooooc;ooooooocooco 

I I I I I t 

CM CM CM CM CM O CM CM CM CM CM CM i-* CM ♦--» CM CM CM 

ooooecooooooocoocc 

CiOOOCQOCDCOOOCOOCOC 
> 4 'LT' 4 '>tcr(X)(Mir' 4 *—‘r^vOirh-«-‘cmf^ 
—i^ao>t^cM— ‘Ov 0 ^in>t»-‘fn> 4 -ca'CT 
mo'rnp^mmP^LnNtcMf^'.o— 
Oso»-'ocoD>tCfvcrcr(Mj.-‘in«-‘irinNtvC 



CC C ocoooo OOOOOOOC c 

I II II I 



m rr CM xp c^ CM m c^v! ro rn m m <n >4- cn > 4 - 
000000000 OOOOOOOC 



cn 



3 

-J 

O 

o 

<r 

CM 



I/) 

JC 

U 

CX 

O 



COCCOOOCQ 

oeoscoocoo 

CMCM^' 0 > 4 "COinO 

>rmooir>roa'inNCLr. 

m ^ ro p-i CM o cr CM 



OCGQQOOO 
or ooocoocM 
CMsOOOOOOCMOin 
ino*'inoofM »-40 
CMsCCM> 4 -fM^{MCM 



000000000000000000 

III I 

rnr <’ i > 3 -> j - c <'( MO (< Trnrornrnfnrr > fn > 4 " 

O OOOOOOOOOCOOOOOO 

O OOOOaOOOOOOOOOQO 
O nTOOCTC coo 00 OOO CMC 00 
— ‘ NTO^n-oC'OocMNTCMvr ooin.cM 
p^ > 4 'sOP^^p^inin.cocroin.o> 4 "'-^o.r- 
— 'Osoo.— i-MP- 4 in<— «^> 4 'CGLr.— i> 4 'nO'— ' 



ooocj>ooooc oooooooooo 

I I I I I III 

p- 

o rn o o o ro o rn o r<^, CM o ro rn po 

CCOOOOCCC v^OOOOOCCO 



o o 
o o 
o o 

o o 

• • 

o o 



CM CM •-* CM CM n J CM CM CM ^ CM CM CM CM (M (M ^ CM 

oc~)Ooeooooooocooooo 

I 

Ln OOC 00 COO OOOOC OOOOC 

^ (Mpw>3-rr'.P^rriOCMinco>ta'r“»^sOvC rnCM 

5: sOCMOOf^rn>rcOvO'n>t cnsOOCMCMCMOLn 

O OO ' CCOCM .- Ma ' ' Oc ^' nc ^ Of <' lODcMp -« cn ^ 4 ' 
«J st >t oo^-^Ln—iin— •in^NCfn' 4 "cn> 4 - oor- 

C 

o 0000 00000 00 000000c 
III I 

• 4 “ 

CM CM O CM CM m CM CM P -4 CM CM CM CM CM CM CM CM CM 
000000000 00 CO OOCOOO 
UO 

jS o qoq 000000 ooo ooo oq 

O ajoocrcMO'oomNr cTcM^p^-a cr^inm^ 
cc nC— ^ r^cMCMCMmcMOOvTinr-ooino^co 
O C ^ P^^ C O > 4 " O c<“‘ <3 O P ^ P ^ cn O o O ' 
o cmcmp^cm— ♦ omo'cn'-^rnvC'O'— ‘incMst^ 



OOOOOOOOOO OOOOOOOC 

I I I I I III 

«k 

> 4 " CM CM CM CM CM CM CM CM CM (M CM CM CM O CM CM CM 

00 00 OOCOOO OOOOOOOC 



ooaooQQOQOoaooooao 

OsO' 4 'Of\lCMOOOOOOOCOOOO 
ocln^ 4 ' 0 >T^omnoxcMooocM 0 ^ 4 '^^m 
cMo^o^x 4 -<^r^'-' 0 <j'>r> 4 '^> 4 -sONOx 
Xr ^ N 4 - cD > 4 - sO ^^ CM > 4 '' 4 - CM > 4 ' cnsOXinp ^ 
CX 

UJ oooooc 00000 00 00000 

cu I II I I I I I 

X 

3 fncn^rn' 4 'rnrororc'rnrororo(Mfn' 4 -' 4 *> 4 ' 

-2: 000000000000000000 

X oocjQQoooooaaoQOQoo 
•- oooop^nC ooooc>ooovOOom 
X moCOCMP^^-OLrUT'XXCMOOrc.CMOCM 
p- h-x insor^x^a. xi 0 WOC 7 o cninxcn 
<L XCM> 4 '*“-‘^OCMCM'n*-‘r-i-‘CM.-ir<^-H^.-i 
X ••••••••••••••••«• 

000000000000000000 

II III I I 



II 


II 


uo 












< 


< 


c 


X 


0 


0 




LU 


z 




-J 


X 


<l 


c 


0 


2: 


-J 




U- 


3 


OJ 


X 


3 




z 






X 






X 


—4 






P-H 


X 


X 


X 


X 


H" 


a 


a 


p— 


< 


x 


X 


< 


X 


X 


X 


2 . 




lU 


LU 







OOOaOQOQOOOOOOOOOO 
xr-o^' 4 'OvOO> 4 'CMa'p^'Tn“^>t-^^o>r 
in X cr vO ^ CM X cM O O' CM CT^ o cm <3 
so o ^o X X in CM o ^ in m ^ o cr 
xc 7 ' XxpcMr <' CMro > 4 ' ro > j - cMn . mxmcMCM 



0000 OOOOOC 00 00 00 00 

I II I I I I I 

mCMO’-^CMCMCMCMCMCMCMCMCM—'CMCMCMCM 

000000000.000000000 

OOOaOQOOOOQQQOOaQCJ 
CM n 4“ X P^ >r o in X. ^<M o CT' m cn cr X >T in 
un > c . ccMP ^.- i *-‘ CMmcMmr ^ r ^ inN 4 “' m — ^ 
o ^ J ^^ X ' 0 (;^ xrno •' Ln ^ 4 '^ ccMX - n ^ c'•ox 
•-* CM in X o in cn in X ^ sO ^ CM — * >t in vO >4" 



000000000000000000 

II III I I 



95 



rn ro rn CO cn (\i rn Nt rn 

oooocooooooooooooo 

ooooocoooooooocooo 
cr 00 sO 0 ^ >4* <M O' iT\ o 'O r<^ ^ o CO o 

CO vO >t o '4' f^J 00 o o f\j h- cr <\j ^ m 

str^ooar^in*-<r^io>tu^a' 0 ^fvjf\j.-ih- 



oooooooooooooooooo 

II I III 

-H cn ro cn cn o ^sJ fsj ro cn fo rr^ CO cn >t cn cn 

oc oc oocooooooooooo 

QOCCCiOQOQOOOOQCCOC 
00 lo r- >t m O' >t fvj O' m ' 4 " fo 

CT' sO <N IT f\j CO o r-4 o 00 >4* o u> 
^ir>^stcr-o»-'vor- 0 '-^o<)NO<Nstr*- 
>Of\jfnrf'roLnO'oo»-H> 4 'vor^*-<»- 4 stc\jr^ o' 



oooooooooooooooooo 

II I I II 

cororncT'>tNtOrr\rvjcnrr'rroof\J 4 'rn> 4 ' 

OOOOOOOOOOOOOOOOOO 

OCOOCOOOQQOOOOOOOO 

^f'-f\jsCr^O'ooO'' 4 “OoaONCOrnojaor- 

CO 00 ^ f\j (\j sT) r- o ro 00 o <Nj o 

ro o o r- rn r- r- CM in nO iT\ <N Nt o r*- ^ 



OOOOOOOOOOOOOOOOOO 
I I I I I I 

rnrni^vt > 4 *'^ rncncHNtO. CM stM- 

OOOOOOOOOOOOOOOOOO 

OOOQQOOCQCOOOOOCOO 
O' h- fMof^r^cMLnvCoinooQoO'str- (M 
O— <r^rno<MM' 0 (MLnocoocoorvj> 4 'f*^ 
m 00 cn 00 rn o (M (M in o 00 r- o n 4" sC >t ^ vO 
OnO*— •»- i'-JOCcn^N 4 “'“<fVf\!N^CO> 4 *»— 



oooooooooooooooooo 

I II II I 

rncOf\J> 4 'rnN 4 'f\jrnfMOrncororocns 4 ‘r<^' 4 ' 

oooooooooooooooooo 



tn- CCCjCQOCQQOOCCQOQQC 

^ O' ^cM' 4 ’CMnO'“<inir.fMs 0 ocr*-if\jv 0 ^ 

3 =. Nrr^'-<r^oaDfMa'vO(Ma'sO*-‘NomvC 2 o>n 

ZD >cco<M > 0 '-^ooNroOfM^inr^r-f\joor 

-j r' > 4 ''^'^fM— <O'(Mr^a 0 fMh*cMcnrvj— <<M^ 

u 

O OOOOOOOOOOOOOOOOOO (M 

III I o 

nT 

fM f^rsjrnrncnor^cnr'orn >T o 

OOOOOOOOOOOC oooooo O 

to o 

QOOOOOOaQQOOCOCOOO o 

LJ m sO sor^M'O'oou'— ^LnLncMvOQOCMOM'fM o 

oc COO-^r-iM-OOOOOvOmcr rXJsOOC^vOO o 

in rn in. vC O' oo<?n 4 t^ o 

vt) — <oooro..— »^^> 4 '»— i^s^'r**-' 4 ’»^ccunf--<^ m. vO 

• •••••••••••••#••• cr CO 

OOOOOOOOOOOOOOOOOO 

I I I I I III 




in 



to 

u cc: 

-J UJ 
O CQ 

o 3 : 

LL o 
> 

X 

X •— > 

— Q 1 

K 

♦— <1 
< >_ 
3 . 



o.r<'rnrornrOf<',ornrc',fnf\jrno^orno 

OOOOCOOOOOOOOOOOOO 

QQOQQQOQQQQOOOQQQO 

ro O' -r in ^ 00 >t <M Nt r- r- CT' M- — • o o 

^^ina''OOinorns 0 ^^roooo'QO<MC 7 ' 

s 0 '“<c'f^Jf^crfMO'O'^fMin> 4 -rnr^r^r*^^ 

xoo^ 4 ’^-' 4 •^ 0 •-ll-••-HM'lnr^JvrmvO^-LAcx 3 



OOOOOOOOOOOOOOOOOO 
I II I I I I I 

corn»-<rO'TM’rnrnrocnfnrnr<'i(Mr^>^> 4 '>^ 

OOOOOOOOOOOOOOOOOO 

QQOaOQQQOQQaOOQQQCJ 
>oovr(no'^u^ii'>x^of'*xiLncMLnvTr^u> 
>OinmTLn^r-sO'HX»^NOrnr>.-<inm^' 
^ O' m m CO o o Lo rn tn r*- O' »-• fM in 00 (M 
0 'fM^'-^.-i^<M(Mm^h-^(M^ro^^^ 



OOOOOOOOOOOOOOOOOO 
II III I I 



II 

to 

2 : 

o 



u 



II 


II 


< 

X 


OJ 


UU 


LU 


2: 


3: 


h- 


N«rt 


►— t 


1— • 




►— 


LL 


a 


o 


o 


LU 


LU 




to 


to 


oc 


LL 


Q. 


LU 


< 


< 


CD 


-J 


-1 


X 


LU 


LU 


O 



96 



OOOO^CVJ 

oooooo 

I 

OQOOOO 

Oc30Ofn^^- 

r^cMc^roro%o 

• ••••• 

oooooo 
I I I 









00 O*^*^ 

OOOOOO 






















1555D 

1642D 

9420D- 

2030D 

75750 

33440- 






















• ••••• 
















t/) 






OOOOOO 
















CL 






1 1 
















o 






















QL 






OOOOviOO 
















QC 






OOOOOO 
















LU 




oo 


oooooo 




in 


ro 


>r 


LT 






UJ 




Z\ 


OrvJOCO-^ 




0 


0 


0 


0 






CL 




5: 


ocmoh*rom 
















< 




3 


(\jfnoo«--<oro 




Q 


Q 


Q 


C 






ZD 




-J 






(T' 


C\i 


rr^ 


0 






o 




O 


• «•••• 




>r 


<r 


cr 


OsJ 






or 




0 


OOOOOO 




rvj 




in 


Ovl 












1 1 1 




CO 


0 


cr 


OQ 






O 




nO 






0^ 


r- 


m 


00 












^O^O^O 




o> 


m 


1—4 


r* 












<300000 




• 


• 


• 


• 




m 






1/1 


1 1 




<r 


c\j 




f\* 




cu 


< 






oaoQOo 














• 


h- 




a 


voosi^cT'om 














rri 


LLi 




or 


'-^rnmo'fvo 
















CD 






r-vOvOcx)>t^^ 




00 


t/) 


uo 












nO 


i>r (NiCT^fn 




3 


-J 


3 


3 


CT 




•f 






• #•••• 




< 


< 


< 


< 


cr 










oooooo 




3 


3 


3 


3 


• 




> 






1 1 




Of 


0 


0 


Of 


0 








•- 






LU 


uu 


LU 


lU 






OO 






OfVjOOOO 
















LLi 






oooooo 






C\i 




<■ 






^sl 




















II 


»-< 






OQQQOO 
















2 _ 






O' rrj (\i (\i (\j 00 


0 












uu 








aoO'rnin>rQO 


0 












2 - 




oO 




ror^ ornxO<\J 




QC 


QC 


QC 


QC 




►—* 


►— I 






,-H r -4 f— • fO ^ r-H 


Q 


LU 


OJ 


LU 


UJ 




1 - 


2: 




QC 




0 


CO 


CD 


CD 


CD 


II 






X 


LU 


OOOOOO 


0 


z: 


5: 


5: 


X 




0 


•» 




CJ 


1 1 1 


0 


0 


3 


3 


3 


LU 


OJ 


vT 


QC 


2: 




0 


2: 


e 


2: 


2 


2 


cn 


• 




3 


-^ 0^000 


0 












Q. 


0 


<3 


-i 


oooooo 


cc 


2: 


2: 




2 




<i 


2: 


a: 




1 


• 


0 


(3 


u 


0 




3 






X 


OQOQQO 


0 


»— 1 


l»4 






O 


3 


»— 


0 


H- 1 


>j'cr'Naom>T 


II 


QC 


QC 


a 


QC 


LU 






UJ 


QC 


sOaO^OOLTi^T' 




LU 


UJ 


LU 


UJ 


t/> 


3 


Cl 


cx 


h- 


(Niroi^OsJUMSk 


<I 




H- 


H- 


h- 


□u 


< 


>— 


l>*4 


<t 




►- 


i-H 






»-4 


C 


h- 


3 


00 


5. 


• #•••• 


UJ 


QC 


QC 


QC 


C3C 


3 


0 


0 


LU 




oooooo 


m 


0 


0 


0 


0 


LU 


h- 




Q 




V \ 1 




1 













97 



\PPSNDIX D 



BRIEF TREATMENT OF RELATED PROBLEMS 
D-1 Case wh^re there are no Restrictions on C, 

In the development of this thesis, starting with the matrix equa- 
tion CU » V, matrix C was restricted to be symmetric and positive definite. 
Suppose now that these restrictions should be removed, and, indeed, let 
it only be required that C have the correct dimensions to relate matrices 
U and V which may have a different number of rows. 

Let us write equation (4) in a modified form 

oC = -Dz. [(cu-v)^z (cu-v) ] (D-1) 

where Z is arbitrary and might be thought of as representing a set of 
weighting factors. If Z is the unit matrix, equation (D-1) reduces to 
equation (4) . 

From equation (D-1) 

oC-ti [U‘^C"2 CU - U‘^C'^2 V-V'^ZCU -►V’^ZV] 

T 

Note that, in general, C C since the requirement of symmetry 
is not made here. 

Using * to denote partial differentiation with respect to any element 
X, as before, one gets 

«*=&[u''c’'*2CU +uVhc*u -uV*ZV-v'’iC*0] 

-th [(cuu'^-vu'^) (c'^*)(z )] (D-2) 
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Thus o(, is made stationary with respect to each element of C if C 



satisfies the equation 

T T 
CUU = VU 

If U is of rank m, where m is the number of rows in U, then 

T T -1 
C = VU (UU ) 

It is interesting to note, that the result is independent of the 
arbitrary weighting matrix Z. 



D-2 Case where C is Orthogonal. 

This problem arises in the treatment of the orientation of a rigid 

body. 



If C is orthogonal, then 



T -1 
C = C = Y 



,-l 



where the notation Y = C is used for convenience in manipulation. 
As before let x be any element of C, and let * denote 



ax • 

Then the relation CY = I, leads to C* = -CY*C, and equation (4) can be 



written as 



so that 



(u'^u-u'^Yv-v'^cu -»-vV) 

(-u‘^Y*v v'^c ) 

= ta[( Y*) (CUV'^C 



(D-3) 



Thus, C< is stationary with respect to variation in any of the ele- 
ments of C, if C satisfies the equation 

T T 

CUV C = VU 

which will be so if C is any of the several values 

c = (vu'^uv'^)^(uv'^)"^ 
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T T 

VU UV is symmetric matrix, and its square root can be obtained as 
described in Section 4, 

These results were obtained by Professor Brock and are given here 
so as to collect in one place information on this general matrix problem* 



100 



INITIAL DISTRIBUTION LIST 



No. 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library 

Naval Postgraduate School 
Monterey, Calif. 93940 

3. Professor John E. Brock 
Mechanical Engineering Department 
Naval Postgraduate School 
Monterey, Calif. 93940 

4. Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, Calif. 93940 

5. Mr. Carl P. Wilson 
Code 452 

General Dynamics/Electric Boat Division 
Eastern Point Road 
Groton, Conn. 03640 

6. Mr. Leslie Kaldor 

Navy Marine Engineering Laboratory 
Annapolis, Maryland 21402 

7. LT Cesar F. Villaran Tapia, Peruvian Navy 
Moreyra y Riglos 769, Pueblo Libre 

Lima, Peru. 

S.A. 

8. Commander 

Naval Ship Systems Command 
Navy Department 
Washington, D. C. 20360 



Copies 

20 

2 

3 

1 

3 



3 

4 

1 



101 



Security Classification 



DOCUMENT CONTROL DATA • R&D 

(Steueity ctmtmilicmtion of body of and indexing annotation mu9t be entered whan tha overetl report ie clmeeitied) 



1. ORIGINATING ACTIVITY (Corpofte author) 

Naval Poatgro^luate Schoo'' 
Monterey, California 93940 



2a. REPORT SECURITY C L ASStFtC A TION 

Unc 



lb. OROUP 



3. REPORT TITLE 

Dotc^nnlnation of OpLimol CompTiaacr^ and Stifx'nenji Matiricen frora nxpGrirneT>tal 
Data 



4* DESCRIPTIVE NOTES (Type of report and incluaive datea) 

Non 



$• AUTHORfS^ (Laet name, firet name, inittaO 

VIU,.‘ ‘XV. 'I s Cn. ’Ot. c. , LT. , Puruviaa Navy 



6. REPORT DATE 

3 . mber IS 67 



7a. TOTAL NO. OF 

102 



7b. NO. OF RKFS 



• a. CONTRACT OR GRANT NO. 



b. PROJECT NO. 



d. 



l i / 



N/A 



$a. ORiatNATOR*S REPORT NUMEER^S; 

1 J / A 



9b. OTHER REPORT NO(S) (Any Other numbera tfia^ may ba aaaijnacf 
Aia tesfort) 

'A/l. 



10. AVA ILABILITY/LIMITATION NOTICES 






11. SUPPLEMENTARY NOTES 



12. SPONSORING MILITARY ACTIVITY 



N/A 



IJavA'al Schcol 



IS. ABSTRACT 



A linear ntructure can bo charactorired by its cuapliance matrix C, which 
in 6x6 cynirnatrical and positive definite and x^hich relates a force 6-vector 
Fj to a displacement 6- vector Dj by the relation CFj - D.. The inverse, iJ, 
of C, is called the stiffness matrix and satisfies Slj This thesis 

deals with the problem of finding optimal values of such matrices C and S 
from^experimental determinations of a sufficient number of vector-pairs 
(Fj,Dj) v/hich are presumed to contain random errors. 

J. E. Brock has introduced this problem area, suggested several different 
criteria of optimality, and solved some of the corresponding spec5,fic pro- 
blems. This thesis completes the solution to a previously unsolved specific 
problem of this group and contributes caiiputationa lly convenient new solu- 
tions to another. Moreover, a computer program, originally written for the 
CDC 1604 has been rewritten, in FORTRAN IV Tianguage, as t\7o program.^ for the 
IBM System 360 computer, and the capability has been significantly augmented 



DD 



FORM 

1 JAN 64 



1473 



IG3 



Unclr-y- ' 

Security Classification 



Unclassified 

'' Security Classification 



I 4. 



link a 



KEY WORDS 



LINK B 



LINK C 



ROLE 



W T 



ROLE 



W T 



ROLE 



W T 



Matrices 
Compliance 
Stiffness 
Analysis of Data 
Experimental Data 
Data Reduction 
Linear Systems 
Least-Squares 
Experiment 
Optimal Matrices 
Matrix Equations 
Matrix Square Roots 



I 




Unclassified 



S/N 0101 -007-6821 



Security Classification 



104 



A-31409 




- jr 




f 



f 



r ‘ 



5 ^ W - 






■-V fh\^^'’"' r i ' 

» 5^ 4 ij|f. -'gi. 



5'?> 



\ 











4 





I 

I 





r 




thesV68 




DUDLEY KNOX LIBRARY 



3 2768 00414859 3 

DUDLEY KNOX LIBRAR'T 





